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Abstract 

A recently introduced effective quantum potential theory is studied 
in a low momentum region of phase space. This low momentum ap- 
proximation is used to show that the new effective quantum potential 
induces a space-dependent mass and a smoothed potential both of 
them constructed from the classical potential. The exact solution of 
the approximated theory in one spatial dimension is found. The con- 
cept of effective transmission and reflection coefficients for effective 
quantum potentials is proposed and discussed in comparison with an 
analogous quantum statistical mixture problem. The results are ap- 
plied to the case of a square barrier. 
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1 Introduction 

In nano-MOSFETs, resonant tunneling diodes and other micro-electronic 
devices, quantum effects are known to play a significant role p. Recently, 
there has been some efforts to create models incorporating quantum effects 
for micro-electronic devices in a less expensive way than, for instance, the 
Wigner-Poisson model [2]. One of these simplified approaches is based on ef- 
fective quantum potentials |3j. The aim of effective quantum potential meth- 
ods is to replace the quantum transport equations by a classical transport 
theory with a modified potential, reproducing to some extent the essential 
quantum effects of the problem. To be faithful, an effective quantum poten- 
tial theory must have similar or lower computational complexity than the 
quantum transport equations. In addition, it must reduce to classical trans- 
port theory for vanishing Planck's constant. Apart from more pragmàtic 
purposes, effective quantum potentials have also attracted attention as an 
attempt to return to phase space and classical statistical mechanics. For a 
general introduction on the subject of effective potentials, the reader can see 
the works [U E] as well as the references therein. 

Recently jïïj|ü], a new effective quantum potential has been proposed and 
applied [7j to model nano-MOSFETs of 25nm. The advantage of the new 
method over the previous ones like is that it does not require the introduction 
of a fitting parameter to be adjusted at the beginning of the simulations. In- 
deed, for effective quantum potentials constructed from Gaussian smoothing 
of the classical potential, there has been some controversy about the correct 
value of the smoothing parameter. Specifically [El E] effective quantum 
potentials Vcaussiç) obtained from Gaussian smoothing of the classical po- 
tential V c (q) according to 



the smoothing parameter being a. The smoothing parameter is related to the 
nonzero size of the electron wave packet and its value changes in a decisive 
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way the position and the value of the peak of charge density in simulations jl] . 
For some candidates for a, see reference jl]. In contrast, in the new effective 
quantum potential method the size of the wave packet is determined by the 
particle's energy, with no need of adjustable parameters. The derivation 
of the effective potential comes from the fact that the Wigner-Bloch and 
quantum Boltzmann equations should have identical stationary solutions. 
The details on the derivation can be found in E] . 

The effective quantum potential introduced in E] depends in a nontriv- 
ial way on momentum. This obscures its interpretation and the construction 
of analytical results. From the conceptual viewpoint, this is a disadvantage 
over methods based on Gaussian smoothing only, where the effective quan- 
tum potential is not momentum dependent. To circumvent these difhculties 
and to obtain a more profound understanding of the new effective quantum 
potential, the present work proposes a low momentum approximation, res- 
tricting to a low momentum region of phase space. As will be seen, for these 
low energy cases the new effective potential amounts to the introduction of a 
space-dependent mass and of a smeared potential. For simplicity, we restrict 
to one-dimensional problems. In this case, the low momentum approxima- 
tion allows for the complete integrability of the equations of motion. From 
the exact solution of the dynamics we can get detailed information about the 
structure of the effective quantum potential theory. 

One of the main objectives of effective quantum potentials is the de- 
scription of tunneling over barriers using a classical language. Neverthe- 
less, tunneling is a nonlocal phenomena which probably can not be exactly 
matched by local effective potential theories 4|. In addition, of course it is 
somewhat misleading to talk about transmission coefficients and so on when 
using effective quantum potentials. Indeed, in this case we would have sim- 
ply "tunneling" as a result of partides undergoing classical transport over 
the barrier. Nevertheless, it is interesting to consider in detail the transport 
of partides in our effective quantum model. For this purpose we consider 
here the specific case of a square barrier, introducing an effective transmis- 
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sion coeficient associated to the effective quantum potential. This effective 
transmission coeficient is analog to the traditional transmission coeficient of 
quantum mechanics, but has a different nature. We define it as the fraction 
of partides surpassing the barrier of the effective quantum potential, for a 
beam of partides going from the left and with equally spaced energies rang- 
ing from zero up to the height of the classical barrier. A similar definition 
is introduced for the effective reflection coeficient. Even if these effective 
coeficients are not strictly analog to the quantum ones, they are useful for 
a quantitative analysis of tunneling in effective quantum potential theories. 

The work is organized as follows. In Section 2, the explicit form of the 
effective quantum potential is introduced. The effective quantum potential 
is expanded considering small energies, leaving us with an approximated 
form characterized by a space-dependent mass and a smeared potential. In 
Section 3, the exact solution for the classical mechanics associated to the 
approximated effective quantum potential is constructed. With the aid of 
the analytical results of Sections 2 and 3, in Section 4 the case of a square 
barrier is analyzed in detail. In particular, we explicitly obtain the effective 
quantum potential and the effective transmission and reflection coeficients. 
We compare these coeficients with the transmission and reflection coeficients 
found in the strictly quantum case. Section 5 is reserved to the conclusions. 



2 Classical Mechanics with Effective Quan- 
tum Potentials 

Our effective quantum potential theory is described by 

• - d£Cff ■ d£&íí (i) 
^ dp ' ^ dq 

where the effective Hamiltonian is 

e eff (q, P ) = ^ + V^(q,p). (2) 
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In equation (J2J), the effective quantum potential is not simply a Gaussian 
smoothing, being defined according to the recent works [SJE|, 

V°"{q,p) = J T(q-q',p)V c (q')dq', (3) 

where V c (q) is the classical potential and 

r(g - q ,p) = - J — smh(— ) exp[-^— + -q)]dk. (4) 

In the last equation, (3 = (k^T) -1 as usual. Observe that the integral 
smoothing (JHJ) reduce fluctuations, a welcomed result especially in the neigh- 
borhood of hetero-junctions. 

From (jïïj) it follows that the effective quantum potential is momentum 
dependent. This fact renders the equations of motion difficult to solve in 
general. However, the form of the integral (J3J) suggests a useful approxima- 
tion. Indeed, let us restrict to small momenta, 

(3hpk/m<t:l, (5) 

to expand the hyperbolic sine at (jU). Of course, the approximation (JSJ has 
to be carefully justified since k is not fixed in the integral, so that, even for 
small p, the condition © may be not satisfied. However, the exponential in 
(JU) is negligible for large k, provi ded 

(3h 2 k 2 /m > 1 . (6) 

When the argument of the hyperbolic sine at is not negligible, so that 
Phpk/m w 1, equation (jüj) can be used to disregard higher order terms 
provided 

/5p 2 /m<l. (7) 

Indeed, substituting Tik w m/((3p) into (jüj), we get ((7|). For even large 
wave numbers the approximation becomes still more accurate. In view of the 
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above arguments, for not too large momenta it is justifiable to expand the 
hyperbolic sine at (jlj) up to 0(f3p 2 /m). This gives 

^,p)^„) + i(^--L)^ + o((^!) 2 ), (8) 

where 

V(q) = dkdq'exv[-^+ik(q-q')]V c (q') 7 (9) 

Í7T J 8m 

M- 1 (q)=m- 1 + dkdq l k 2 exp[-^—+tk(q-q')]V c (q') . (10) 

24:%m 2 J 8m 

Notice that the expansion is to all orders of h, the only restriction being to 
consider small momenta. Taken into account (|ÏÏJ), we obtain the approximated 
effective quantum Hamiltonian 

etlU = ^ + V(q), (11) 

the Hamiltonian for a particle of position-dependent mass under a time- 
independent potential. Quantum corrections are present both in V(q) and 
M(q). 



3 Exact Solution 

Hamilton equations for £ e JJ prox (q 1 'p) rea d 



p . dV 1 dM 2 

9= m' p = -ir q + 2M 2 i[ q - p ■ (12) 

Eliminating p and restricting to level sets of constant energy e e JJ prox = e, 
the result is 

.. _ l_dV_ (V - e) dM 

q ~~M~dq~ + M 2 ~dq~' { } 

or 

m = —£+ , (14) 
6 



with the potential 

Tïl 

V Q (q) = J ^(V(q)-e)+e. (15) 

The interpretation of ()14j) is as follows. At least for low momenta, the effec- 
tive quantum potential induces modificat ions of the original classical poten- 
tial, so that the final result is classical motion under the potential (|15p. The 
exact solution for ()14|) then follows from elementary methods. For H — > 0, we 
have M — > m and V — > V c , so that V® — > V 7 and the classical Newtonian 
equation is recovered. The addition of the irrelevant numerical constant e 
at the end of (|Tïï|) was just a matter of convenience to obtain the correct 
classical limit when fi — ► 0. 

The usefulness of the form (|15p together with (|Ü HTÜ|) is that it provides a 
clean way to observe the modifications produced by quantum effects in the 
effective quantum potential model. For a given classical potential, we have 
a recipe for constructing the potential V® and to study the corresponding 
completely integrable Newtonian equation. In particular, this strategy can 
be used to obtain conclusions about tunneling rates related to the smearing 
of the classical potential when replaced by the potential V Q . This possibility 
will be explored in the next Section. 

4 Tunneling 

Consider a square potential barrier centered at q — 0, height Vq and width 
2L, given by 

V c = V [9(L + q) + 6{L -q)-l], (16) 

where 9 is the unit step function. The potential (fTïï|) was chosen for its 
simplicity only; the specific choice of the form of the potential barrier does 
not affect in a decisive way the results. Using (p iTÜ|) . we obtain 
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where erf is the error function. The quantum potential V® trien follows 
from (j!5jl . 

Graphics of the space-dependent mass are shown in figure 1, with m — 1, 
L = 1/2, Vo = 1, = 1/8 and two different vàlues of h, namely h = 10 (full 
line) and H = 30 (dotted line), using arbitrary units. As expected, M — > m 
for |ç| 3> L. However, in the neighborhood of the barrier there is a small 
increase of the mass. In contrast, inside the barrier the mass is lowered. For 
larger quantum effects, we observe a displacement of the location of the mass 
peak, as well as an increase of the mass peak. 

In figure 2, we consider the potential V® for m — 1, L — 1/2, Vq = 1, 
f3 — 1/8, £ — .25 and for varying h. For h — 0, the classical square barrier 
is recovered. For larger quantum effects (h = 3 for the dashed line and 
h = 6 for the dotted line), we see an increasing smoothing of the potential 
V® . The height of V® becomes smaller for larger quantum effects, increasing 
tunneling. The observed smoothing is somewhat reminiscent of the smearing 
of the self-consistent hole potential for quantum plasmas described by the 
Wigner-Poisson system [TOj- In spite of the small increase of the mass M(q) 
near the barrier, the net result of the effective quantum potential is clearly 
the lowering of the barrier. 

To obtain a condition for tunneling, consider that, for not too big quan- 
tum effects, the maximum of V® is obtained for q = 0. For very large 
quantum effects, we can show graphically that the maximum of V® is not 
at q = anymore, being displaced to symmetric positions around the ori- 
gin. However, in these cases the smoothing is very strong. Disregarding this 
possibility, we obtain the analytical expression 
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x \V o erf\^0^-e\+e (19) 



For a classical particle of energy e to surpass the barrier V®, it is necessary 
that e > max(V Q ). With (fTÏÏj) . this gives 

* V . .20) 



the tunneling condition in our classical description. For — > it reduces sim- 
ply to e > Vq as expected. Also, tunneling is increased for small temperatures 
and for a small thickness of the barrier. 

As said before, in the present classical model for a quantum process the 
word "tunneling" is to be used with caution. However, we can even assign 
an effective "transmission coemcient" to be calculated with the smeared out 
potential V®. Indeed, let a beam of classical partides going from the left 
against the barrier V c . Suppose that the beam is constituted by partides 
with energies ranging from to Vo, with a homogeneous distribution in en- 
ergy. In other words, the probability to find a particle with energy in a range 
5 E is simply 5 E/Vq. Classically, no one of the partides of this hypothetical 
beam would surpass the barrier. However, if we admit the model of the effec- 
tive quantum potential (jïïj) and restrict to small energies, so that j3V <C 1, 
we can apply the approximation of the preceding sections and define an ef- 
fective transmission coefficient t given by the ratio of the number of partides 
that surpass the barrier to those that are reflected. In view of (j2(Jj) . this 
effective transmission coemcient is 

V.-Voerf(^) N2mL 



Similarly, we obtain an effective reflection coemcient given by 

r = erf ( . (22) 
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Of course we have t + r = 1. Graphics of t and r are shown in figure 3, where 
once again we see that larger quantum effects increase tunneling. In this 
particular example, quantum effects are measured by the non dimensional 



quantity H = J p/m h/L. Significant tunneling begins to happen for H m 2. 
For H = 2.97 we have already 50% of the partides surpassing the barrier. 

We have defined a transmission coefhcient well suited for the classical 
description we have adopted. It is interesting to compare this with the 
truly quantum description. The quantum version of the preceding situa- 
tion would be a statistical mixture of incoming eigen states of well defined 
energy, {ip(x) = e tkx }, where ip(x) is a particular eigen function on the en- 
semble. Here, we are considering partides going from the left (k > 0) with an 
energy E = h 2 k 2 /(2m). Hence, we have dE/Vo = h 2 k dkj (mVo), so that the 
probability of finding a state with wave-number k in the ensemble would be 
p(k) = h 2 ko/ (mVo), for < k < k and zero otherwise, where k = y/2mVo/h. 
We have the proper normalization J fc ° p(k) dk = 1 as it should be. 

What is the "transmission coefficient" associated to this statistical mix- 
ture? For just one incoming eigen state e ífcx , any text book on quantum 
mechanics gives the result 

t Q = A: 2 (A; 2 - fcg) 

° ne k 4 - k 2 k 2 - {k 4 /4) sinh 2 (2 v /A; 2 - k 2 L) 

for this transmission coefficient. For the statistical mixture just described, 
due to the linearity of the Schròdinger equation the transmission coefficient 
is given by the average 

t Q = / t one (k) p{k) dk , (24) 
J o 

or, after some easy scalings and rearrangements, 

t Q = 2 í 1 x\x 2 -l)dx 

Jo x 4 -x 2 -(l /4) sinh 2 (2 yi^VQ) ' 1 ; 

In the quantum transmission coefficient is a function of the parameter 
Q = [kç^L)^ 1 = h/ (y/mVo L) only, while the parameter H = ^J/3/mh/L is 
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present in the effective quantum potential description. In figure 4 we show 
the quantum transmission coeficient v* as well as the quantum reflection 
coefíicient = 1 — v* as functions of Q. For Q = 1.75 we have already 
= 0.5. As seen comparing figures 3 and 4, the transmission coefficients for 
the effective quantum potential and the quantum statistical mixture problem 
have similar functional behaviors in spite of the somewhat different analytical 
expressions (|2Tj) and ()25j) and the different relevant quantum parameters, H 
and Q. In a sense, this supports the viewpoint that V Q adequately simulates 
tunneling, even in our low momentum approximation. Other examples (dou- 
ble square barriers, other classical potential forms) can be easily constructed, 
with results similar to those of this Section. 

5 Conclusion 

In this work, the effective quantum potential (jïïj) was studied in a low energy 
region of phase space, characterized by /3p 2 /m <C 1. In its original form, the 
effective quantum potential is a complicated function of momentum. In the 
low momentum approximation, we identified a space-dependent mass and a 
smoothed potential, both obtainable from the classical potential. For one- 
dimensional problems, the dynamics of a particle with space-dependent mass 
under a time-independent potential is completely integrable, a fact that al- 
lows for a number of conclusions. In particular, we have found the potential 
V® at equation í)15|l. which, in the low momentum approximation, describes 
the classical motion under the effective quantum potential. In connection 
to the problem of a square barrier, we proposed effective transmission and 
reflection coefficients for effective quantum potentials, which were explicitly 
found with the aid of the exact solution for the equation of motion. Compar- 
ison with true quantum transmission and reflection coefficients is allowed if 
an appropriated quantum statistical mixture is introduced. In the case of a 
square barrier, we identified two non dimensional parameters measuring the 
relevance of quantum effects for tunneling, H = J (3/m Ti/L for the effective 
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quantum potential model and Q = %/ (y/mVo L) for the quantum statistical 
mixture. 

In conclusion, we have obtained a more detailed physical intèrpret ation 
of the effective quantum potential (jïïj), in terms of a space-dependent mass 
and a smoothed potential, in a low momentum approximation. It would 
be interesting to generalize the results of this work to more dimensions. In 
particular, it would be interesting to investigate the integrability of potentials 
V® obtained from higher dimensional integrable classical potentials. 
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Figure 1: Space-dependent mass with m — 1, L — 1/2, Vo = l,/3 = l/8 and 
two different vàlues of h, namely H = 10 (full line) and H = 30 (dotted line), 
using arbitrary units. 



Figure 2: The potential V Q for m = 1, L = 1/2, Vb = 1, /3 = 1/8, e = .25 
and ^ = (full line), = 3 (dashed line) and % = 6 (dotted line). We 
consider arbitrary units. 



Figure 3: Transmission coeficient t (full line) and reflection coefficient r 
(doted line) as functions of the non dimensional parameter H = 



Figure 4: Transmission coefficient t Q (full line) and reflection coefficient r Q 
(doted line) as functions of the non dimensional parameter Q = h/(\fmV ü L). 
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